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Spin qubits defined in carbon nanotube quan- 
tum dots are of considerable interest for encod- 
ing and manipulating quantum information be- 
cause of the long electron spin coherence times 
expected. However, before carbon nanotubes can 
find applications in quantum information process- 
ing schemes, we need to understand and control 
the coupling between individual electron spins 
and the interaction between the electron spins 
and their environment. Here we make use of spin 
selection rules to directly measure - and demon- 
strate control of - the singlet-triplet exchange 
coupling between two carbon nanotube quantum 
dots. We furthermore elucidate the effects of 
spin-orbit interaction on the electron transitions 
and investigate the interaction of the quantum 
dot system with a single impurity spin - the ulti- 
mate limit in spin sensitivity. 

A powerful method to probe the spin dynamics of 
quantum dots is by measuring electron transport in a 
double quantum dot device in the spin blockade regime 
In this transport regime, the tunneling of an electron 
between the two quantum dots is forbidden by spin selec- 
tion rules and hence the current is suppressed. However, 
spin blockade can be lifted by the interaction of the elec- 
tron spins with their environment and a measurement of 
the (leakage) current thus directly probes these interac- 
tions [2-4]. The main spin relaxation and decoherence 
modes in carbon nanotube that have been considered so 
far are hyperfine and spin-orbit coupling. The impor- 
tance of hyperfine interaction in the electron spin dynam- 
ics of carbon nanotubes has recently been demonstrated 
in 13 C enriched nanotube quantum dots [5|,l6|. The sig- 
nificance of spin-orbit coupling, a result of the curvature 
of the carbon sheet, was demonstrated for carbon nano- 
tubes [THIOj and is expected to be relevant for graphene 
as well [9(. 

A further important consideration in any realistic de- 
vice is the presence of impurities or defects and their 
coupling to the electron spins. For example, magnetic 
catalyst particles are used in nearly all carbon nanotube 
synthesis methods while defects such as vacancies are 
also widely regarded _to giv e rise to magnetic behavior 
in carbon materials 
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While defect densities in 
high-quality carbon nanotubes are extremely low, their 



presence is still highly relevant in the context of quantum 
information processing in which the coupling to even a 
single impurity spin will affect device performance. A 
single impurity spin, however, is extremely difficult to 
detect. Here we achieve this by using the exceptional 
sensitivity of the transport current of a double quantum 
dot in the spin blockade regime to the spin environment. 
We show that the interplay of an impurity spin and spin- 
orbit interaction has a dramatic effect on the spin states 
of the double dot and find excellent agreement with a 
theoretical model. 

The device we consider is a single- walled carbon nano- 
tube grown by chemical vapor deposition using methane 
with natural isotope ratios and contacted by Au elec- 
trodes. Side and top barrier gates are used to define and 
control the double quantum dot, see Fig. 1(a). A typ- 
ical charge stability diagram of the device is shown in 
Figs. l(b,c) in which the ordered pairs (n,m) indicate 
the effective electron occupancies of the many-electron 
double quantum dot. In the presence of a source-drain 
bias voltage of V s d = — 1 mV, a honeycomb structure, 
characteristic of a double quantum dot [TUTU is clearly 
visible. The large-small-large-small alternation of the ad- 
dition energy in the honeycomb pattern indicates that the 
electron states are spin degenerate but that the orbital 
degeneracy of the nanotube has been broken. 

Of particular interest is the observation of Pauli spin 
blockade |2(j, of which a characteristic feature is the 
strong bias dependence of the current for every other 
added electron as seen in the top and bottom rows of 
Figs. l(b,c), respectively. This is further illustrated by 
the detailed measurements in Figs. l(d,e) that corre- 
spond to the region in Figs. l(b,c), highlighted by the 
dashed circle, in which the excess number of electrons 
changes from (0,1) to (1,2). For negative bias voltage 
the current is strongly suppressed while a large current 
is observed at the base of the bias triangles when the bias 
is positive. 

These measurements can be understood considering 
that, for negative bias, a flow of electrons from the left 
to the right quantum dot necessarily involves a transition 
from the (1,1) to the (0,2) charge state. Since the (0,2) 
ground state has to be a singlet by virtue of the Pauli ex- 
clusion principle, the (l,l)-»(0,2) transition is forbidden 
by spin selection rules when the electrons on the dou- 
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Figure 1: (a) Schematic of the carbon nanotube double quantum dot device used in the experiment, (b) Charge stability 
diagram of the double dot for V s d = — 1 mV. The circles indicate the regions of spin blockade, (c) Charge stability diagram of 
the same region for V B d = +1 mV. (d,e) Detailed measurements of the bias triangles at the (0,1) to (1,2) charge transition for 
both polarities. Note the applied bias of 2 mV which allows us to observe the first excited states in the bias triangles, and the 
difference in current scales for the two figures. Since the inter-dot capacitance is small, the bias triangles for the electron and 
hole cycle strongly overlap. The energy level diagrams illustrate possible and forbidden transitions in the double quantum dot. 



ble dot form a T(l,l) triplet state, see also schematics 
in Figs. l(d,e). When the bias is positive the S(0,2) — »• 
S(l,l) transition involves singlet states only and electrons 
can freely move from the right to the left quantum dot 
and no current suppression is observed. 

In the following we make use of spin selection rules to 
directly probe the spin system of the double quantum dot 
in detail. Our main results are illustrated in Fig. 2. The 
top row shows the stability diagrams of spin blockaded 
bias triangles for five different barrier potentials. The 
bottom row shows the current through the double dot as 
a function of magnetic field B and detuning e for three 
fixed barrier potentials, Vbar =0,10, and 15 mV, corre- 
sponding to three of the stability diagrams in Fig. 2(a). 
Note that in Fig. 2(b) we plot the normalized current to 
accentuate the evolution of the spin states by effectively 
subtracting the background current. The complete mea- 
surement set, for all measured barrier voltages, as well 
as examples of measurements at other bias triangles are 
shown in the Supplementary Information. 

As evident from the measurements presented in 
Fig. 2(b), the observed evolution of the spin states can be 
very rich and, as discussed below, deviates considerably 
from naive expectations based on the simple even-odd 
spin filling pattern observed in Fig. 1. We start with a de- 



scription of the current dependence at Vt, ar = mV, cor- 
responding to the leftmost stability diagram in Fig. 2(a). 
In the Pauli blockade regime only the (1,1) triplet states 
are occupied and no current can pass through the de- 
vice. However, when the magnetic field is non-zero the 
degenerate triplet splits and for a given magnetic field 
and detuning the triplet energies equal the hybridized 
singlet energies, see schematics in Fig. 3(a). In the pres- 
ence of spin relaxation such as due to spin-orbit medi- 
ated electron-phonon interaction [To| . the spin flip rate 
is much higher at these points and provides an efficient 
escape route from the blockaded triplet state. Thus mea- 
suring the current as a function of detuning and magnetic 
field maps out the exact energy dependence of the sin- 
glet energies, and therefore the singlet-triplet exchange 
energy J(e), as seen in the leftmost funnel-shaped pat- 
tern in Fig. 2(b). We obtain an excellent fit to the ex- 
perimental data using a theoretical model of the device 
(for details see Methods section) which yields a tunnel 
coupling t — 70 /ieV, consistent with an independent es- 
timate of t from the stability diagram in Fig. 2(a). The 
measurements also demonstrate our ability to electrically 
tune J(e) by varying the detuning which, in the limit of 
large negative detuning, approaches J(e) ~ 4t 2 /U, where 
U ^10 — 20 meV is the charging energy of the quantum 
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Figure 2: (a) Bias triangles at the (0,1) to (1,2) charge transition for V s d = —1 niV and B = T for five different barrier gate 
voltages, (b) Normalized current as a function of detuning and magnetic field for Vb ar = 0, -10, and -15 mV. The detuning 
(energy difference between the (1,1) and (0,2) charge states) axes follow the dashed white lines in panel (a) for the three 
respective barrier gate voltages. In the leftmost panel two sets of curves of high current are visible. In the middle panel an 
additional set of curves appears close to zero magnetic field. The rightmost panel shows a series of additional faintly visible 
current peaks as illustrated by the linetrace and the inset which shows part of the measurement (\B\ < 0.37 T) with enhanced 
contrast. The evolution with detuning and magnetic field of the features closely follows that predicted by the model described 
in the main text as illustrated by the dashed curves. The free fitting parameters are the interdot tunnel coupling and the 
strength of the Heisenberg interaction between the double dot and a spin-1/2. 



dots, see schematics in Fig. 3(a). The calculated current 
for this situation is shown in the leftmost plot of Fig. 3(c). 

While the funnel-shaped pattern in the leftmost panel 
in Fig. 2(b) has been observed previously in other sys- 
tems such as GaAs the evolution of the spin states in 
the remaining two panels has not previously been seen in 
any double quantum dot. When the barrier voltage is set 
to Vbar =10 mV a completely different pattern appears in 
which two sets of curves that approach a common asymp- 
tote at B ~ 0.26 T are seen, see dashed-dotted line. In- 
triguingly, when the barrier voltage is increased further 
to Vbar =15 mV, the two sets of curves are accompanied 
by three further, weaker, sets of curves [see the rightmost 
panel in Fig. EJ^b)]. 

To explain our observations, we propose a model where 
additionally to the electrons on the double dot a third 
spin, the impurity spin, is present. As justified below, 
we assume it to be a spin-1/2 that couples to the spin of 
one of the dot electrons via an isotropic Heisenberg inter- 
action with a coupling strength Ji mp , see Fig. 3(b). The 
strength of the interaction as well as that of the tunnel 
coupling are tuneable by the gate electrodes. The rele- 
vant states of the combined quantum dot and impurity 
system can be characterized by their total spin: A four- 
fold degenerate spin 3 /2 state and two doubly degenerate 
spin 1/2 states, as shown in Fig. 3(b). The S = 3/2 states 



cannot mix with the S — 1/2 (0,2) state since tunneling 
conserves spin, and therefore block the current, while all 
the spin 1/2 states can take part in transport through 
the device, having a (0,2) component. The multiplets 
split in a magnetic field with the energy of states with 
higher magnetic spin quantum numbers mj passing the 
energy of states with lower magnetic spin quantum num- 
bers as indicated in Fig. 3(b). At finite magnetic fields, 
the lowest lying state in the 5 = 3/2 quartet (Q4) has the 
highest occupation probability. Since selection rules limit 
the possible transitions to those from state Q4 to D2 and 
£>4, we therefore expect to see two strong curves, tracing 
out the shape of the 5=1/2 levels. This is indeed the 
case as observed in the data, see Fig. 2(b), which is in 
excellent agreement with the calculated current as shown 
in Fig.3(c). 

The model also allows us to investigate the effects of 
spin-orbit interaction on the electron transitions. A first 
indication of the presence of spin-orbit coupling is the 
zero-field dip in the spin blockade leakage current around 
e = as seen in the leftmost plot in Fig. 2(b). This fea- 
ture has previously been observed in carbon nanotube 
[f| and InAs [2l[ double quantum dots and has been ten- 
tatively attributed to spin-orbit interaction. These re- 
sults were reproduced in recent theoretical work in which 
spin-orbit interaction was shown to introduce non spin- 
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Figure 3: (a) Energy of the relevant two-electron states of a double quantum dot as a function of detuning. The black line 
represents the (1,1) triplet states (S = 1), while the red and orange curves represent mixtures of a (0,2) singlet state with a (1,1) 
singlet state (5 = 0). For non-zero magnetic field the triplet state energies split. For a given detuning and magnetic field there 
are points where the triplet energies equal the singlet energies (marked by dots), (b) The energy of the relevant two-electron 
states of a double quantum dot in the Pauli blockade regime including a spin- 1/2 impurity as a function of detuning. The 
impurity spin couples via isotropic Heisenberg interaction with a coupling strength Jimp to one of the dot spins, see inset. 
The black line represents a S = 3/2 quartet. These states are blocked and will have high occupation. The remaining curves 
are S = 1/2 doublets. Spin selection rules allow transitions with Am a = ±1: Q3 -o- Dl, Q3 «-» D3, Q4 «-» D2, Q4 o D4 
and D2 D3, while other, forbidden, transitions are Q4 Dl and Q4 <-> D3. (c) Current as a function of magnetic field 
and detuning. Left: for no impurity spin (Jimp = peV) and t = 70 /ieV, Middle: for impurity spin with Jimp = 8 ^ieV, and 
t = 87/^eV, no spin-flip tunneling, Right: for impurity spin with Jimp = 6/ieV, and t — 170 peV, in the presence of spin-flip 
tunneling with t s = 7.5 fieV. Note that the origin of the curve D2 «-> D3 in the rightmost plot is the large interdot tunnel 
coupling and not the spin-orbit interaction. 



conserving tunneling between the two quantum dots [22j . 
In our model this is characterized by a spin-flip tunneling 
amplitude t s . Since in the presence of spin-orbit interac- 
tion m s is no longer a good quantum number, spin selec- 
tion rules can be violated, resulting in additional tran- 
sitions and therefore extra curves in the current plots. 
These curves are indeed observed in the data, as most 
clearly seen in the rightmost plot in Fig. 2(b). The posi- 
tion and evolution of the additional three curves are in ex- 
cellent agreement with those predicted by the model cal- 
culation as illustrated by the dotted curves in Fig. 2(b) 
and the corresponding calculated current in Fig. 3(c). 

Spin-orbit interaction also affects the shape of the 
curves. The spin-flip tunneling induces coherent tran- 
sitions between T + /_(l,l) and S(0,2) resulting in an 
avoided crossing between the triplet and singlet states, 
the size of which depends on the (0,2) component of the 
singlet. This results in relatively narrow inner and a 
wide outermost curve which increases in width as e — > 0. 
A comparison with the data gives fairly narrow bounds 
for the spin-flip tunneling and yields t s — 7.5 /ieV, see 
Supplementary Information. Using the estimate from 
Ref. [22] and an orbital energy E ur i, of the order of 2 



meV we can deduce the spin-orbit interaction energy as 
Eso — E or bt s /t w 0.1 meV. We note that this value 
is similar to previous estimates @, 0] even though the 
strength of the spin-orbit interaction is not a priori clear 
in the many-electron limit and for mixed orbital states. 

The above analysis demonstrates that we are able to 
detect the presence of a single impurity spin coupled to 
the carbon nanotube double quantum dot and determine 
its spin quantum number. An important question that 
remains is the nature of the impurity. Even a single mag- 
netic atom absorbed on the nanotube wall would have 
a significant effect on the spin states [23| and could be 
present as a result of the nanotube growth process, see 
Methods section. Further possibilities are charge traps 
in the gate oxide or the presence of a defect such as 
a vacancy or dopant in the carbon lattice. Even high- 
quality carbon nanotubes have been found to contain 
one defect per 4 /im on average [24[. This would im- 
ply that the majority of our devices (having total length 
L ~ 1/im) contain either zero or one impurity. Whereas 
many nanotube devices will therefore be without a single 
defect, requirements on the nanotube quality will become 
more stringent when many quantum dots are coupled in 
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a large-scale quantum circuit. 

Our experiments also show that the exchange inter- 
action between the carbon nanotube and impurity spin 
can be precisely controlled with a gate electrode, see Fig. 
2(b). In combination with the recent advances in attach- 
ing single atoms or molecules to the nanotube sidewall 
|25|, this suggests the possibility of storing quantum in- 
formation into the attached groups and using the carbon 
nanotube as a quantum bus and for spin state read-out 
[26l . 27 1 . The ability to control these interactions will be 



instrumental in developing carbon materials for quantum 
information processing. 



Methods 
Device fabrication 

Carbon nanotubes were synthesized by chemical va- 
por deposition using a procedure similar to Ref. [28| . To 
form the catalyst particles, ~ 0.5 /xg/cm 3 Fe(N03).9fi20 
(Sigma-Aldrich) was sonicated in isopropanol. The bare 
Si/Si02 substrate (300 nm thermal oxide) was dipped 
into the solution immediately after sonication and dried 
by air blowing. It was then heated to 900 °C in a 1" 
tube furnace under a hydrogen gas flow of 400 seem dur- 
ing which the iron nitrate was reduced to iron particles. 
Once at this temperature, 500 seem of methane (con- 
taining natural isotope ratios) was added to the flow for 
8 minutes. After the methane flow was stopped, the sub- 
strate was cooled under the same 400 seem hydrogen flow. 
Alignment marks were defined by electron-beam lithogra- 
phy and scanning electron microscopy and atomic force 
microscopy were used to locate the nanotubes in rela- 
tion to these. Further steps of electron-beam lithogra- 
phy and evaporation contact the nanotubes with 20 nm 
thick, 300 nm wide gold ohmic contacts, separated by 700 
nm, which also form the outer barriers of the quantum 
dots. The final electron-beam lithography layer writes 
the 100 nm wide and 100 nm spaced gates, which are 
made from two 1.2 nm layers of aluminum, oxidized in 
air and capped with a further 20 nm of titanium and 6 
nm of gold. The source lead is grounded by the virtual 
earth of a xlO 7 current to voltage preamplifier and the 
drain has a small variable dc bias applied. The chip is 
mounted on the cold finger of a 60 mK dilution refriger- 
ator and connected to the measurement apparatus with 
filtered low-frequency lines. 



Theoretical model 

We use a master equation technique to model the phys- 
ical system of interest that consists of a double quantum 
dot coupled to an impurity spin. We are interested in the 



electrical current through the double quantum dot in the 
presence of spin relaxation and decoherence. 

The dynamics of the electrons on the double quantum 
dot is modeled by an extended two-site Hubbard Hamil- 
tonian with single orbitals. The impurity spin is modeled 
as a spin- 1/2 coupled to one of the dots via an isotropic 
Heisenberg interaction. To facilitate the numerical eval- 
uations we represent the density matrix of the system in 
a many-body basis resulting in 32 basis states. The cal- 
culations are performed in a rotated basis in which the 
Hamiltonian of the system is diagonal. 

The double quantum dot is connected to left and right 
metallic leads which are described as reservoirs of non- 
interacting fermions. The applied magnetic field does not 
affect the distribution function of the fermions, as long 
as the Zeeman energy is small compared to the Fermi 
energy. The tunneling between the leads and the double 
quantum dot is treated perturbatively to second order in 
the tunneling amplitude. To model spin relaxation and 
decoherence each spin is coupled to an independent bath 
of harmonic oscillators in thermal equilibrium with a flat 
spectral density. The coupling between the spin and the 
bath is treated perturbatively to lowest non-vanishing 
order in the coupling strength. 

Within the Born and Markov approximations we derive 
an equation of motion for the density matrix of the double 
quantum dot and the impurity spin, that we solve for 
the stationary state. Using a similar methodology we 
derive the current operator up to second order in the 
tunneling amplitude and calculate the electrical current 
in the stationary state. Further details of the model can 
be found in the Supplementary Information. 
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Experimental data 

In the following sections we present additional data to that in the main text. Figures Sl(a,b) show stability diagrams 
of the device used to obtain the various capacitance terms and charging energies of the two quantum dots following 
the work of van der Wiel et al []| . The parameters are summarized in Table I. The difference in the charging energies 
U, as well as the single-particle energies AE might be related to geometry: QD# is slightly smaller than QDl, see 
inset in Fig. SI (a) which shows an SEM image of the device taken after measurement. 

Figures Sl(c,d) show measurements for various magnetic fields (B = 0,4,8 T) perpendicular to the nanotube axis 
at both negative and positive applied bias. These measurements provide additional evidence for the interpretation of 
the data in terms of spin blockade and allow us to identify excited states in the bias window. Fig. Sl(c) shows that 
at finite field the current at the lower edges of the spin blockade triangles are suppressed as expected when the triplet 
states split off and spin exchange with the leads is suppressed. The spin-blockade region also becomes smaller as the 
T_ states move down in energy with field and transitions to T(0, 2), indicated by the arrows, are observed at smaller 
detuning. As shown in Fig. Sl(d), the behavior is precisely the opposite for positive (V s d = +2 mV) bias voltages. At 
B = T the double dot is not spin blockaded and a large current at the base of the triangles is observed. However, 
as the field increases, the double dot enters the spin blockade regime as the T_ triplet pushes the base of the triangle 
out, but blocks the S(0, 2) T-(l, 1) transition, see schematic in Fig. Sl(e). 

Figure S2 illustrates how the normalized data as presented in the main text is obtained from the current as- 
measured. Fig. S2(a) shows the current as a function of detuning and magnetic field for spin-blockaded bias triangles 



Table I: Capacitance terms and charging energies of the double quantum dot 



Property 


QDl 


QDr 


Unit 


Comment 


C 


15 


6.2 


aF 


Quantum dot capacitance 




0.22 




aF 


Interdot capacitance 




0.12 


0.060 


aF 


Plunger gate-QD capacitance 


u 


11 


26 


meV 


Charging energy 


U' 


0.37 




meV 


Electrostatic coupling energy 


AE 


1.8 


2.0 


meV 


Single-particle addition energy 
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Figure S 1: (a) Charge stability diagram for V B d = —1 mV. The arrows indicate the spacing in Vuftlright °f the hexagonal 
pattern. Inset: Scanning electron micrograph of the device. The arrow indicates the position of the carbon nanotube. (b) 
Detailed measurement of the bias triangles at the (0,1) to (1,2) charge transition (effective electron occupancies) for V B a = —2 
mV. (c) Bias triangles for V sc i = —2 mV and B = 0, 4, and 8 T. The arrow in the measurements at B — 4, 8 T indicates 
transitions to the T(0, 2) state, (d) Bias triangles for V B d = +2 mV and B = 0, 4, and 8 T. Spin blockade is observed at 
large magnetic fields, (e) Energy schematics of the double quantum dot. The symbols correspond to the positions in the bias 
triangles in panels (c,d). 

at V s d = —1 mV and Vb ar = 10 mV, see Fig. 2 in main text. As shown in Fig. S2(c), the current varies strongly 
as a function of detuning. The corresponding normalized current, shown in Fig. S2(b), accentuates the evolution of 
the spin states by effectively subtracting the background current which allows for a detailed comparison with theory 
(dashed red lines). The values on the detuning axis can be obtained from the stability diagrams, as well as from fits 
to the theory. Both methods give very similar estimates for the detuning. In the measurements presented here, the 
values on the detuning axes are those obtained from fits to the theory. 

We have observed a dependence of the current measurements on the direction and rate of the magnetic field sweeps. 
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Figure S 2: (a) Current as-measured as a function of detuning (projected on Vi e ft) and magnetic field, (b) Normalized current 
corresponding to the section of panel (a) between the dashed lines. The dashed red curves are fits to the theory, see main text, 
(c) Current as function of detuning (projected on Vi e ft) at B = T. (d) Current as function of magnetic field for Vi e ft = —1.0 
V. (e) Current as function of magnetic field for Vuft = —1.1 V for both sweep directions, as indicated by the arrows. 



This is most clearly seen for positive detuning, as illustrated in Fig. S2(e) for a sweep rate of 20 T/h. Note, however, 
that the position of the current peaks in magnetic field appears to be unaffected. As a nearby Ge-thermometer showed 
similar behavior we attribute the observed hysteresis to an artefact of the measurement system (possibly heating) and 
not to coupling of the electrons to the nuclear spin system (our nanotubes contain ~ 1% 13 C). 

Figure S3 shows the complete normalized data set, measured at V s( i = — 1 mV for nine different Vbar, for the 
spin-blockaded bias triangles labeled A in Fig. Sl(a). The coupling to the impurity spin is small in Fig. S3(a) but is 
clearly observed in Figs. S3(b,c). For relatively large t, additional curves (dashed blue lines) become faintly visible, 
see Fig. S3(c). All measurements allow for a detailed fit to the theory with the results for t and Ji mp summarized in 
Fig. S4. Both the tunnel coupling between the quantum dots and the coupling to the impurity spin strongly depend 
on - that is, are tuneable by - the barrier gate electrode. The tunnel coupling increases for increasing Vbar, most 
strongly so after Vbar = 12.5 mV while the coupling to the impurity spin increases and then decreases again. 

Finally, in Fig. S5 we show measurements on two sets of spin-blockaded bias triangles different from those discussed 
in the main text. The stability diagrams in Fig. S5(a,b), both measured for Vbar = mV, correspond to the regions 
labeled B and C in Fig. SI (a), respectively. The results are very similar to those discussed in the main text, 
demonstrating that our results are reproducible and that, apart from the parity of the electron occupation (even or 
odd), the qualitative features of the measurements do not depend on the total number of electrons on the quantum 
dots. We furthermore verified that the observed features and the evolution of the spin states with detuning and 
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Figure S 3: (a-c) Measurements of the normalized current as a function of magnetic field and detuning for nine different Vbar- 
In panel (a) the exchange interaction J; mp is small and only a single set of lines is visible. In panel (b) the exchange interaction 
is stronger and an additional set of lines appears. The evolution of the lines can be fitted very well with a theoretical model 
of the device, see dashed red lines. For the measurements shown in panel (c), the tunnel coupling between the two quantum 
dots is relatively strong and three additional, faintly visible, sets of lines appear, as indicated by the arrows. As described in 
the main text, the evolution of these lines follows that predicted by the theoretical model, see dashed blue lines. 

magnetic field do not depend on the applied bias (tested for V s d = — 0.5, —1, —2, —3, —5 mV). As expected, the 
features are not observed for non-spin blockaded bias triangles. Also note that in Fig. S5(a), we observe a strong 
coupling of the electrons on the quantum dots to the impurity spin {Ji mp = 8.5 /ieV), even for Vbar = mV. This 
indicates that the coupling to the impurity spin is tuneable not only by the barrier gate electrode (as shown in the 
main text), but also by the side gates. 
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Figure S 4: Evolution of the quantum dot-impurity coupling strength Ji mp (left) and tunnel coupling t (right) with barrier gate 
voltage V bar - 
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Figure S 5: (a,b) Spin-blockaded bias triangles and spin-state evolution as a function of detuning and magnetic field for two 
spin configurations different from those discussed in the main text. The dashed red lines are fits to the theory. 
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Theoretical Model and results 



In this section we describe in detail the theoretical model and present results for the electrical current for some 
special cases. The total system consists of a double quantum dot, a nearby impurity spin-1/2, the metallic leads, and 
bosonic heat baths. This system is modelled by the Hamiltonian 

Htot = H.S + Hleads + Ht + Hb + HsB- (1) 

The system Hamiltonian, H$, models the double quantum dot and the impurity spin, Hi eai j s (Hb) models the leads 
(heat baths) and Ht (Hsb) models the interaction between the leads (heat baths) and the system. Specifically, we 
have 

Hs = H DD + H M + Hi, (2) 

where Hdd is a Hubbard Hamiltonian describing the double quantum dot, Hm is due to the applied magnetic field 
and Hj models the interaction of the double quantum dot with the impurity spin. For the double quantum dot we 
have 

2 2 

H D d = ^2e l n l - y~] W(ci g c 2g ' + c\ a ,c 1<y ) + Ujn^nn + U'nm 2 , (3) 

i— 1 <7.<r f i—1 

that allows up to two electrons per dot. The number operator is rii = X)a nicr = Scr C L C ^ f° r dot i — {1,2} and 
spin a = {t,4-}- The operator c\ a (ci a ) creates (annihilates) an electron on dot i with on-site energy Si. The spin 
conserving tunnel coupling between the two dots is t aa > = t for a = a', and the spin- flip tunnel coupling between the 
two dots is t aa i = t s for a ^ a' ' . Ui is the charging energy of each dot and U' is the interdot Coulomb energy between 
the dots. The Hamiltonian part due to the applied magnetic field is 

i=0 

where i — refers to the impurity spin and the spin operators are defined in the standard way er^ = J2 aa , c^a-^^da' , 
where a is the vector of the 2x2 Pauli matrices. Aj = g^sB is the Zeeman splitting due to the magnetic field B 
along z, g- factor gt and the Bohr magneton [is- We consider a Heisenberg interaction between dot 1 (i = 1) and the 
nearby impurity spin (i = 0) of the general form 



Hi = JxO-qO-i + J y o\o\ + J z GqUI, 



(5) 
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with J x , J y , J z being the strength of the interaction, and we set J x = J y = J z = J imp for isotropic interaction. 

Having described the system of interest, i.e., the double quantum dot and the impurity spin- 1/2, we turn to the 
metallic leads, the bosonic heat baths and their interaction with the system. The left and right leads are described 
by a Hamiltonian of the form 

Hleads = "^2 e ikd eka dika, (6) 
ika 

where d\ ka (deka) creates (annihilates) an electron in lead I = {L, R} with momentum k, spin a and energy eek- The 
interaction between the dots and the leads is given by the tunneling Hamiltonian 

R t = ^2(tLc\ a d LkrT + t R cl a d Rka ) + H.c, (7) 

ka 

where t^tn) is the tunnel coupling between dot 1(2) and lead L(R) and we consider the symmetric case where 
tL = tR. To include spin relaxation and decoherence we consider a generic bosonic heat bath that is modelled as a 
set of harmonic oscillators and is described by the Hamiltonian 

2 

^B = EE ^i,j a h a i,3- ( 8 ) 
i=0 j 

We assume that the impurity spin and each dot are coupled to an independent bosonic heat bath (i — refers to 
the impurity spin bath and i = 1, 2 to dot 1, 2) and there are no environment-induced correlations between them. 
Moreover, we assume that the spins are flipped one at a time. The operator a\ - (ciij) creates (annihilates) a boson 
in mode j and Uij are the frequencies of the bath modes. The impurity spin and the spins on the double quantum 
dot interact with the corresponding bath via the general model Hamiltonian 

2 

• H - ( -- ( 9 ) 

i=0 j 

where the spin-flip operators are = cj, and A^ j is the coupling constant between the impurity spin (i = 0), 
dot 1 (i = 1) or dot 2 (i = 2) and the jth mode of the corresponding bath. Hsb allows spin-flip processes for the 
impurity spin and the spins on the double quantum dot, via energy exchange with the heat bath, which lead to a 
leakage current in the spin blockade regime. 

To investigate the electronic transport through the double quantum dot we employ a master equation approach 
and derive an equation of motion for the reduced density matrix, p, for the system of interest (double quantum dot 
and impurity spin). Starting with the density matrix of the total system, Xtot, and within the standard Born and 
Markov approximations we derive an equation of motion for p by tracing over the leads and bosonic baths degrees of 
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freedom, i.e., p=TiE{xtot} where Tr^j...} means trace over the environmental degrees of freedom. It can be shown 
that the equation of motion of the density matrix p can be written in a compact form with the superoperators £g, 
Cieads and L B as 



p{t) = C S p{t) + Cleadsp{t) + C B p(t), (10) 



with the free evolution term 



C sP (t) = ~[H s ,p(t)}, 



and the terms due to the electronic leads 



1 f°° 

Ci eads p{t) = -^Tr B { / dr[H Tl [U{r)H T W \r), p(t) ® Pleads}}}, 



and the bosonic baths 



1 f°° 

C B p{t) = -^Tr £ { J dr[H S B, [V{t)H SB V\t), p{t) ® p B }}}. 

Both Ht and Hsb are treated in second order perturbation theory with the operators U (r) = exp[— i(Hs+Hi eac [ s )T/h} 
and V(t) — exp[— i(H$ + Hb)t/H}, and pi ea ds, Pb being the equilibrium density matrix for the leads and the bosonic 
baths respectively which have the same temperature. 

The quantity of interest is the electrical current that flows through the double quantum dot. By definition the 
current operator, for example, for the right lead is 

Ir = -eN R = -e % -\H T ,N R } = e l -Y^{tR(^ !T dR k a-t R S Rka c 2a ), (11) 

ka 

where Nr = ^2 ka d Rka duka is the number operator for electrons in the right lead. Tracing out the leads we derive 
the expectation value of the current, (Ir), that is calculated in the stationary state, i.e., when p = in Eq. (|10p . in the 
sequential tunneling regime. The left lead current is {Il)=-(Ir}- For the double quantum dot with single orbital levels 
there are sixteen possible states and the maximum number of electrons is four. Including the two possible impurity 
spin states results in a total of thirty two basis states. For the calculations we write Eq. tfTU]) as a set of coupled 
equations for all the matrix elements of p, in the energy basis, including the normalisation condition for the diagonal 
elements that express the occupation probabilities. The resulting system, for p = 0, is then solved numerically. The 
theoretical model predicts for the current a rather small increase in the range |Aj| < ksT/A, due to spin relaxation, 
and this can be understood through the Boltzmann factor that arises naturally within the model. 
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Figure S 6: As Fig. 3(c) (middle) in the main article but for different (/-factors for the impurity spin and g\ — g^ — 2 for the two 
dots. 



13 5 7 



13 5 7 



g -0.5- 




Figure S 7: As Fig. 3(c) (middle) in the main article but for anisotropic Heisenberg interaction between the impurity spin and 
one of the quantum dots. 



Finally, we examine the electrical current as a function of magnetic field and energy detuning for some interesting 
cases. Figure £j6] shows the current when the impurity spin and the double dot spins have different g-factors. As the 
difference increases the curves of high current begin to split and this effect becomes noticeable at about go = 0.8(/i, 
with <?i = <?2 = 2 for the two dots. The exact pattern of the curves can be identified through the dependence of 
the three-spin energy spectrum on magnetic field and detuning when the double dot is in the spin blockade regime 
(see also Fig. 3(b) in the main article). In order to reproduce the experimental data with the theoretical model the 
difference in the (/-factors should be less than about 20%. 

Figure £[7]shows the current for anisotropic Heisenberg interaction between the impurity spin and one of the dots. We 
consider antiferromagnetic coupling (J x , J y , J z > 0), which produces a pattern for the high current curves consistent 
with the experimental data. Ferromagnetic coupling (J x , J y , J z < 0) produces, in general, a different pattern. As 
the anisotropy increases additional curves appear, similar to the case of spin-flip interdot tunneling described in the 
main article (see Fig. 3(c)). However, the analysis of the functional dependence of the curves on magnetic field and 
detuning shows that it is impossible to fit the theoretical results of high anisotropy to the experimental data. Our 
calculations suggest that the anisotropy has to be much smaller than 20%. 

Figure £[5] shows the current for different spin-flip tunnel couplings between the two quantum dots. As explained 
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in the main article with increasing coupling additional high current curves appear. Further, the width of the curves 
increases with the spin-flip coupling. This is attributed to coherent transitions between the three-spin states, which 
are induced solely when the spin-flip coupling is nonzero. These transitions act in addition to incoherent transitions 
induced by spin-flips because of the coupling to the heat baths. Both effects define a limit on the magnitude of 
the spin-flip coupling observed in the experiments. For a too weak coupling no additional curves would be observed, 
whereas for a too strong coupling the curves would be washed out as shown in Fig. £[5] The theoretical model indicates 
that the spin-flip coupling should be on the order of t s ~ 10 /xeV and the best fit to the experimental data is obtained 
for t s w 7.5 fj,eV. 
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